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ABSTRACT 

Store  separation  testing  is  undertaken  to  ensure  that  stores  can  be  safely  released  from  aircraft 
in  flight.  Grid  testing  in  wind  tunnels  produces  data  that  can  be  used  to  predict  store 
separation  behaviour.  This  document  describes  the  algorithms  used  to  prepare  test  plans  for 
store  separation  grid  testing  and  to  reduce  grid  data  to  the  form  required  for  analysis.  The 
information  is  presented  in  a  general  form,  and  may  also  provide  a  useful  reference  for  other 
forms  of  wind  tunnel  testing. 
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Executive  Summary 


The  release  of  stores  from  an  aircraft,  including  fuel  tanks  and  ordnance,  poses  a  risk 
that  must  be  managed  in  order  to  safely  employ  stores  during  operations.  Under  some 
conditions,  particularly  at  high-speed,  the  store  may  make  unintended  contact  with  the 
aircraft,  endangering  the  lives  of  personnel  and  causing  equipment  damage. 
Management  of  the  risk  involves  ground  and  flight  testing  to  understand  the 
behaviour  of  stores  during  separation  from  the  aircraft.  Wind  tunnel  testing  is  an 
integral  part  of  stores  clearance  projects. 

The  Defence  Science  and  Technology  Organisation  (DSTO)  has  an  integral  role  in 
supporting  Royal  Australian  Air  Force  (RAAF)  stores  clearance  campaigns  through 
wind  tunnel  testing  and  store  trajectory  prediction.  Store  separation  tests  for  modern 
fighter  aircraft  are  undertaken  in  the  DSTO  transonic  wind  tunnel.  DSTO  uses  a  grid 
technique,  where  aerodynamic  data  are  gathered  in  a  prescribed  grid  pattern  below  the 
store  carriage  position.  This  grid  covers  the  region  that  the  store  is  expected  to  travel 
through  after  release  from  its  parent  aircraft.  The  aerodynamic  data  is  then  combined 
with  other  information,  including  mass  and  moments  of  inertia  of  the  store,  to  predict 
the  store's  trajectory  following  release  from  the  aircraft. 

The  algorithms  required  to  prepare  a  test  plan  for  a  transonic  wind  tunnel  grid  test  are 
described  in  this  document.  Also  included  are  the  processes  used  to  reduce  wind 
tunnel  grid  data  to  the  form  required  for  analysis  of  store  release  trajectories.  The 
implementation  of  these  algorithms  for  specific  applications  is  not  discussed,  but 
examples  are  provided  to  aid  with  development  and  validation  of  software.  The 
algorithms  are  presented  in  a  general  form  and  this  document  may  provide  a  useful 
reference  for  other  forms  of  wind  tunnel  testing. 
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Abbreviations  and  Acronyms 

AIAA  American  Institute  of  Aeronautics  and  Astronautics 

AOSG  Air  Operations  Support  Group 

ASC  Aircraft-Stores  Compatibility 

BMC  Balance  Moment  Centre 

c.g.  Centre  of  gravity 

DSTO  Defence  Science  and  Technology  Organisation 

DSTOres  Defence  Science  and  Technology  Organisation  release  evaluation  suite 
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MATLAB  MATrix  LABoratory 

MRC  Model  Moment  Reference  Centre 

RAAF  Royal  Australian  Air  Force 

SS  Store  Support 

STEME  Store  Trajectory  Estimation  in  a  MATLAB  Environment 

TWT  Transonic  Wind  Tunnel 


Angles 


(f>c,b/  0c,  b /  y/c,b 

(f>ac,c/  0ac,c ,  yfac,c 

$g,ac,  @g,ac/  yfg,ac 
(/)g,b,  0g,b,  ySg,b 


$g,ss/  Og,ssr  y/g,ss 

&  &  £ 


a,  P 


Roll,  pitch  and  yaw  angles  of  the  store  body  axes  relative  to  carriage 
axes. 

Roll,  pitch  and  yaw  angles  of  the  carriage  axes  relative  to  aircraft 
axes. 

Roll,  pitch  and  yaw  angles  of  the  aircraft  axes  relative  to  tunnel  axes. 
Roll,  pitch  and  yaw  angles  of  the  store  body  axes  relative  to  tunnel 
axes. 

Roll,  pitch  and  yaw  angles  calculated  from  feedback  of  the  store 
support  arm  sensors. 

Angular  misalignment  between  balance  and  store  body  axes,  defined 
in  Appendix  B.l. 

Incidence  (tangent  definition)  and  sideslip  (sine  definition) 
orientation  angles,  relative  to  the  wind  direction.  In  the  absence  of 
flow  angularity  corrections  this  corresponds  to  the  tunnel  gravity 
axes. 


Displacements 


Xc,b/  yC/b,  Zc,b 
Xgts  })g,br  Zg,b 
Xgr  Vgf  Zg 


Xg,0  ))g,c,  Zg/C 
Xg,SS/  })g,ss,  Zg,SS 

Xb/  yb/  Zb 


Displacement  of  the  store  in  carriage  axes. 

Displacement  of  the  store  in  tunnel  axes,  relative  to  carriage  position. 
Displacement  of  the  store  in  tunnel  axes,  relative  to  tunnel  reference 
point. 

Displacement  of  the  carriage  position  in  tunnel  axes. 

Displacement  calculated  from  feedback  of  the  store  support  sensors. 
Displacement  of  the  model  moment  reference  centre  relative  to  the 
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balance  moment  centre  in  balance  axes. 

xcg/  ycg,  zcg  Displacement  of  the  real  store  centre  of  gravity  with  respect  to  the 

model  reference  point,  in  store  body  axes. 


Forces  and  Moments 


Ht  = 


BALT 


M 


BALT . 


XBALT 


F, 


YBALT 


F 


ZBALT 


M 

M 

M 


XBALT 


YBALT 


ZBALT 


H^  = 


Ft 


BALG 


M 


BALG 


F 


XBALG 


Fx 


YBALG 


Fr 


ZBALG 


M 


XBALG 


M 

M 


YBALG 


ZBALG  . 


H  = 


Ft 


BAL 


M 


BAL  . 


F 


XBAL 


YBAL 


Fr 


ZBAL 


M 

M 

M 


XBAL 


YBAL 


ZBAL 


1 _ 

Mx 

F  = 

Ft 

,  M  = 

My 

fz_ 

_MZ_ 

F 

rXc.g. 

Mx,g. 

Fc,  = 

FYc.g. 

>  M,,= 

Myc.g. 

_FZc.g._ 

Mz,g. 

Tare  loads  in  balance  axes. 


Gross  loads  in  balance  axes. 


Aerodynamic  loads  in  balance  axes. 


Aerodynamic  loads  (forces  and  moments)  in  store 
body  axes. 

Aerodynamic  loads  (forces  and  moments)  about 
the  real  store  centre  of  gravity. 


Aerodynamic  Force  and  Moment  Coefficients 

Cx,  Cy,  Cz  Axial,  lateral  and  normal  force  coefficients  in  store  body  axes. 

C/,  Cm,  Cn  Rolling,  pitching  and  yawing  moment  coefficients  in  store  body  axes. 
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Strain  Gauge  Balance  Parameters 


R 

Ri 

Ro 

R'o 

Rt 

R' 

[Cl] 

[C2] 

H 

H* 

kxr  fcy,  kz, 
k-Mxi,  fcvrx2/ 

kMYb  k-MY2/ 
klMZl,  k-MZ2 

K 

ki,  k2,  . . ./  kio 

x',  y ',  z' 
v,  r|,x 


Vector  of  absolute  strain  gauge  balance  output. 
Component  i  of  strain  gauge  balance  output. 

Buoyant  zero  strain  gauge  balance  output. 

Wind-off  zero  strain  gauge  balance  output. 

Tare  load  strain  gauge  balance  output. 

Wind-on  strain  gauge  balance  output. 

Linear  calibration  coefficients. 

Non-linear  calibration  coefficients. 

Vector  of  measured  loads. 

Vector  of  squares  and  cross  products  of  component  loads. 
Tare  weight  coefficients. 


Matrix  of  sting-balance  deflection  coefficients. 
Sting-balance  deflection  coefficients. 

Axial,  lateral  and  normal  deflections  in  balance  axes. 
Roll,  pitch  and  yaw  angular  deflections  in  balance  axes. 


Transformation  Matrices 


Mc,b  Euler  rotation  matrix  from  store  axes  to  carriage  axes. 

Mac,c  Euler  rotation  matrix  from  carriage  axes  to  aircraft  axes. 

M*  ac  Euler  rotation  matrix  from  aircraft  axes  to  tunnel  axes. 

MgrC  Euler  rotation  matrix  from  carriage  axes  to  tunnel  axes 

Mg/b  Euler  rotation  matrix  from  store  axes  to  tunnel  axes. 

Mg,ss  Euler  rotation  matrix  from  store  support  arm  to  tunnel  axes. 

MSS/b  Pseudo-Euler  rotation  matrix  from  store  support  arm  to  store  body 

axes. 

Mb,bai  Non-Euler  rotation  matrix  from  balance  axes  to  store  body  axes. 

MW/b  Euler  rotation  matrix  from  store  body  axes  to  wind  axes. 

MPrb  Euler  rotation  matrix  from  store  body  axes  to  missile  axes. 

Wind  Tunnel  Parameters 


M  Freestream  Mach  number. 

q  Freestream  dynamic  pressure. 

Pt  Tunnel  total  pressure. 

Ps  Tunnel  static  pressure. 

Re  Reynolds  number. 

Y  Ratio  of  specific  heats. 

p  Freestream  density. 

R  Gas  constant. 

Ts  Tunnel  static  temperature. 

p  Freestream  kinematic  viscosity. 
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V  Freestream  velocity  vector 

u,  v,  w  Components  of  freestream  velocity  in  tunnel  axes. 

Model  Parameters 

d  Model  reference  dimension. 

S  Model  reference  area. 


vi 
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1.  Introduction 


When  stores  are  released  from  an  aircraft  they  may  encounter  a  highly  non-uniform  flow 
field  produced  by  the  parent  aircraft,  particularly  at  high  speeds.  Stores  include  any 
removable  object  that  may  be  attached  to  an  aircraft,  and  which  can  be  released  during 
flight.  This  includes  fuel  tanks,  avionic  and  storage  pods,  missiles,  bombs,  pylons  and 
racks.  The  aerodynamic  interaction  between  the  store  and  the  parent  aircraft  may  lead  to 
contact  between  the  store  and  the  aircraft,  endangering  the  lives  of  personnel  and  causing 
severe  equipment  damage.  For  this  reason,  it  is  important  to  predict  the  separation 
behaviour  of  stores  to  ensure  they  can  be  safely  employed  in-flight.  The  Air  Operations 
Support  Group  (AOSG)  is  the  responsible  agency  for  ensuring  that  the  weapon  is  safe  to 
carry  and  release  from  Royal  Australian  Air  Force  (RAAF)  aircraft.  Australia  uses  an 
approach  based  on  the  methodology  of  MIL-HDBK-1763  [1]  for  certification  of  aircraft- 
stores  compatibility  (ASC)  and  flight  clearances. 

Assessment  of  employment  and  jettison  characteristics  for  ASC  is  one  element  of  the  MIL- 
HDBK-1763  methodology.  For  this  assessment,  the  Defence  Science  and  Technology 
Organisation  (DSTO)  and  AOSG  use  a  'grid  method'  [2]  whereby  an  aerodynamic  force 
and  moment  database  is  generated  by  testing  a  model  of  the  store  in  a  wind  tunnel  facility 
with  a  model  of  the  parent  aircraft.  The  store  model,  mounted  on  a  six  degree-of-freedom 
support  arm  via  an  internal  strain  gauge  balance,  is  positioned  at  pre-defined  grid  points 
in  the  region  where  the  store  is  expected  to  travel  when  it  is  released.  Aerodynamic  forces 
and  moments  are  measured  at  each  grid  point  for  a  range  of  store  attitudes.  These  forces 
and  moments  are  then  used  in  a  trajectory  prediction  computer  program  together  with  a 
database  of  the  freestream  aerodynamics  of  the  store  in  isolation  and  a  variety  of  other 
inputs  including  store  mass  and  moments  of  inertia  data,  ejector  force  and  distribution 
models  and  atmospheric  properties.  DSTO  uses  an  in-house  developed  program,  DSTO 
release  evaluation  suite  (DSTOres),  for  trajectory  prediction,  and  AOSG  use  their  own 
prediction  program.  Store  Trajectory  Estimation  in  a  MATLAB  Environment  (STEME)  [3]. 

Wind  tunnel  testing  for  store  release  at  DSTO  is  conducted  in  the  DSTO  Transonic  Wind 
Tunnel  (see  Appendix  A  for  basic  information  on  this  facility).  A  sub-scale  model  of  the 
store  is  mounted  on  a  six  degree-of-freedom  actuated  support.  This  support  system  allows 
the  store  to  be  rotated  at  varying  roll,  pitch  and  yaw  angles  and  translated  axially,  laterally 
and  vertically  within  a  three-dimensional  plane.  A  sub-scale  half-model  of  the  parent 
aircraft  is  mounted  to  a  turntable  support  on  the  tunnel  test  section  sidewall  enabling  a 
pitch  rotation.  The  validity  of  this  method  has  been  verified  in  a  number  of  store 
separation  tests  [4].  The  store  is  positioned  at  pre-defined  grid  points  in  the  region  it  is 
expected  to  travel  through  after  release  from  the  parent  aircraft.  Aerodynamic  forces  and 
moments  are  measured  at  each  grid  point,  for  a  variety  of  store  orientations,  by  a  strain- 
gauge  balance  that  is  integral  to  the  support  arm. 

This  document  describes  the  algorithms  required  to  prepare  transonic  wind  tunnel  grid 
test  plans  and  to  reduce  wind  tunnel  grid  data  to  the  form  required  for  input  into  store 
trajectory  prediction  software.  Sections  2  and  3  provide  the  definitions  of  axis  systems  and 
angles,  respectively,  used  during  aerodynamic  testing  at  DSTO.  Section  4  describes  the 
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transformations  between  coordinate  systems  and  the  process  of  defining  the  store  grid  for 
testing.  Section  5  outlines  the  process  for  reducing  measured  strain  gauge  balance  data 
into  the  aerodynamic  forces  and  moments  required  for  store  release  prediction.  Example 
calculations  for  both  test  preparation  and  data  reduction  are  presented  in  Appendix  C,  and 
can  be  used  to  validate  software  implementations  of  the  data  reduction  algorithms 
presented  in  this  report. 
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2.  Definition  of  Axis  Systems 

The  following  axis  systems  are  used  in  this  document.  These  axis  system  definitions 
generally  conform  to  the  AIAA  Guide  Nomenclature  and  Axis  Systems  for  Aerodynamic  Wind 
Tunnel  Testing  [5]. 

2.1  Tunnel  Gravity  Axis  System 

The  gravity  axis  system  is  an  Earth  fixed  axis  system,  which  has  its  z-direction  aligned 
with  the  gravity  vector.  The  origin  is  located  at  the  nominal  centre  of  the  test  section, 
which  corresponds  to  the  centre  of  the  sidewall  turntable,  and  is  400  mm  laterally  from  the 
inner  sidewall.  It  is  defined  as  station  5257  in  [6,  §2.1. 2.2].  The  axes  are  shown  in  Figure  1 
and  are  defined  as  follows: 

Xg  Gravity  longitudinal  axis,  perpendicular  to  the  gravity  vector  and  contained 

in  a  plane  parallel  to  the  support  system  pitch  plane,  positive  upstream. 

Yg  Gravity  lateral  axis,  perpendicular  to  the  gravity  Xg-Zg  plane,  positive 

direction  determined  by  the  positive  Xg  and  Zg  directions  in  conjunction 
with  the  right-hand  rule. 

Zg  Gravity  vertical  axis,  collinear  with  the  gravity  vector,  positive  toward  the 

tunnel  floor. 


Figure  1  Tunnel  gravity  axis  system 
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2.2  Balance  Axis  System 

The  balance  axis  system  origin  is  fixed  at  the  balance  moment  centre  (BMC)  about  which 
forces  and  moments  were  resolved  during  calibration.  This  is  defined  as  the  geometric 
centre  of  the  gauged  elements  as  shown  in  Figure  2. 

The  Xbai,  Ybai ,  and  Zbai  coordinate  axes  rotate  with  the  strain  gauge  balance  and  are  defined 
as  follows: 

Xbai  Collinear  with  the  strain  gauge  balance  longitudinal  axis,  positive  towards 

the  metric  end. 

Ybai  Perpendicular  to  the  strain  gauge  balance  Xbai-Zbai  plane,  positive  direction 

determined  by  the  positive  Xbai  and  Zbai  directions  in  conjunction  with  the 
right-hand  rule. 

Zbai  Collinear  with  the  gravity  axis  when  the  strain  gauge  balance  is  orientated 

with  zero  pitch  and  roll,  positive  downwards. 


Figure  2  Balance  axis  system  (DSTO-BAL-08  shown) 

2.3  Store  support  axis  system 

The  store  support  axis  system  defines  the  position  and  orientation  measured  by  the  store 
support  system  encoders.  The  origin  of  the  store  support  axis  system  is  located  at  a  fixed 
distance  from  the  store  support  pitch  pivot  point  along  the  store  support  x  axis,  as  shown 
in  Figure  3.  This  distance  represents  the  chosen  centre  of  orientation  for  the  model,  which 
may  or  may  not  be  coincident  with  the  model  reference  centre  or  balance  moment  centre. 
Measurements  in  the  store  support  axis  system  are  uncorrected  for  balance-sting 
deflection.  The  axes  XSS/  Yss  and  Zss  are  defined  as  follows: 
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Xss  Store  support  longitudinal  reference  axis,  positive  upstream. 

Yss  Perpendicular  to  the  store  support  Xss-Zss  plane,  positive  direction 

determined  by  the  positive  Xss  and  Zss  directions  in  conjunction  with  the 
right-hand  rule. 

Zss  Store  support  vertical  reference  axis,  collinear  with  the  gravity  axis  when 

the  store  support  is  orientated  with  zero  pitch  and  roll,  positive 
downwards. 


2.4  Aircraft  Axis  System 

The  aircraft  axis  system  is  a  body  axis  system  which  has  its  origin  at  the  aircraft  reference 
point.  The  aircraft  axis  system  is  shown  in  Figure  4.  Only  the  directions  of  the  axes  are 
required,  and  these  are  defined  as  follows: 

Xac  Aircraft  longitudinal  reference  axis,  positive  toward  the  nose  of  the  aircraft. 

Yac  Aircraft  lateral  reference  axis,  perpendicular  to  the  aircraft  Xac-Zac  plane  and 

positive  as  defined  by  a  right-handed  system. 

Zac  Aircraft  vertical  reference  axis,  parallel  to  and  in  the  same  direction  as  the 

gravity  Zg  axis  when  the  aircraft  is  upright  and  level  in  pitch  and  roll. 
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Aircraft 


Figure  4  Aircraft  axis  system 

2.5  Store  Body  Axis  System 

The  store  body  axis  system  origin  is  designated  as  the  model  moment  reference  centre 
(MRC),  as  shown  in  Figure  5.  The  axes  Xb,  Yb,  and  Zb  are  defined  as  follows: 

Xb  Model  longitudinal  reference  axis,  positive  toward  the  nose  of  model. 

Yb  Model  lateral  reference  axis,  perpendicular  to  the  body  Xb-Zb  plane  and 

positive  as  defined  by  a  right-handed  system. 

Zb  Model  vertical  reference  axis,  parallel  to  and  in  the  same  direction  as  the 

gravity  Zg  axis  when  the  model  is  upright  and  level  in  pitch  and  roll. 


Xb,  X 

P 

Figure  5  Body  axis  system,  in  black,  and  missile  axis  system,  in  red 
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2.6  Wind  Axis  System 

The  origin  of  the  wind  axis  system  is  coincident  with  the  body  axis  system  origin  at  the 
MRC.  The  axes  XW/  Yw  and  Zw  are  defined  as  follows: 

Xw  Wind  longitudinal  axis,  parallel  to  the  total  velocity  vector,  differs  from 

body  Xb  axis  by  the  angles  a  and  ft. 

Yw  Wind  lateral  axis,  perpendicular  to  the  wind  Xw  axis  and  contained  in  the 

Xw-Yb  plane.  Differs  from  the  body  Yb  axis  by  the  angle  (3. 

Zw  Wind  vertical  axis,  perpendicular  to  the  to  the  wind/body  axis  Xw-Yb  plane 

and  contained  in  the  body  axis  Xb-Zb  plane.  Differs  from  the  body  Zb  axis  by 
the  angle  a. 

2.7  Missile  Axis  System 

The  origin  of  the  missile  axis  system  is  coincident  with  the  body  axis  system  origin  at  the 
MRC.  The  missile  axis  system  rotates  with  the  model  through  yaw  and  pitch  rotations 
only.  The  Yp  axis  is  coincident  with  the  Yb  axis  with  the  store  at  zero  roll  angle  (</>cs)-  This 
feature  has  led  to  the  missile  axis  system  being  called  the  non-rolling  body  axis  system. 
The  axes  Xp,  Yp,  and  Zp  are  shown  in  Figure  5,  and  are  defined  as  follows: 

Xp  Missile  longitudinal  axis,  coincident  with  and  in  the  same  direction  as  the 

body  Xb  axis. 

Yp  Missile  lateral  axis,  perpendicular  to  the  missile  axis  Xp-Zp  plane  and 

positive  as  defined  by  a  right-handed  system. 

Zp  Missile  vertical  axis,  perpendicular  to  the  body  Xb  axis  and  contained  in  the 

(missile  axis)  pitch  plane  defined  by  the  body  Xb  axis  and  an  intersecting 
gravity  vector. 

2.8  Carriage  Axis  System 

The  origin  of  the  carriage  axis  system  is  coincident  with  the  store  MRC  when  the  store  is 
positioned  at  carriage,  as  illustrated  in  Figure  6.  The  origin  is  fixed  with  respect  to  the 
parent  aircraft. 

Xc  Carriage  longitudinal  axis,  collinear  with  and  in  the  same  direction  as  the 

model  Xb  axis  when  the  store  is  positioned  at  carriage. 

Yc  Carriage  lateral  reference  axis,  perpendicular  to  the  carriage  Xc-Zc  plane  and 

positive  as  defined  by  a  right-handed  system. 

Zc  Parallel  to  the  carriage  ejector  stroke  direction,  positive  in  the  direction  of 

the  stroke. 
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Positions  within  the  carriage  axis  system  are  often  defined  in  polar  coordinates  using  the 
following  parameters: 


r 

r/D  (Calibre) 

0p 

(Polar  Angle) 


Location  of  the  store  relative  to  carriage  axes  as  a  radial  distance 
from  the  store  release  point  in  the  Yc-Zc  plane,  as  shown  in  Figure  6. 
r  normalised  by  store  diameter,  D. 

Location  of  the  store  relative  to  ejector  axes  as  an  angle  away  from  Zc 
in  the  Yc-Zc  plane.  Positive  towards  positive  Yc  axis,  as  shown  in 
Figure  6.  (Note:  this  is  opposite  to  the  right  hand  rule  about  Xc.) 
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3.  Definition  of  Angles 


The  following  definitions  of  aerodynamic  and  orientation  angles  are  used  for  aerodynamic 
testing  at  DSTO,  including  tests  performed  in  the  transonic  wind  tunnel. 


3.1  Aerodynamic  Angles 


Testing  in  the  DSTO  wind  tunnels  uses  the  tangent  definition  for  incidence  and  the  sine 
definition  for  sideslip.  Incidence  is  defined  as  the  angle  between  the  Xb  axis  and  the 
projection  of  the  wind  vector  on  the  Xb-Zb  plane. 


w 

tan  a  =  — 


(1) 


u 

Sideslip  is  the  angle  between  the  wind  vector,  V,  and  its  projection  on  the  model  Xb-Zb 
plane,  v, 


sin  B  =  T~~r 
\V\ 

These  angles  are  limited  to  the  ranges 

-  n  <  a  <  n 
-n  <  <7t 

3.2  Orientation  Angles 


(2) 

(3) 


The  orientation  of  one  set  of  axes  relative  to  another  is  commonly  represented  by  Euler 
angles.  Testing  in  the  DSTO  wind  tunnels  uses  the  standard  yaw  (i//)  pitch  (0)  and  roll  (<f>) 
angles,  applied  in  the  order  (y/,  0,  </>).  To  avoid  ambiguity,  the  angles  are  limited  to  the 
following  ranges: 

-  n  <  y/  <  n  (4) 

-E<q<?L 
2  2 

-n  <(j)<n 


3.3  Transformation  between  Aerodynamic  and  Orientation  Angles 


The  aerodynamic  angles  describe  the  orientation  of  the  store  body  axes  relative  to  the 
tunnel  flow  axes.  In  the  absence  of  flow  angularity  corrections,  tunnel  flow  axes  are 
coincident  with  tunnel  gravity  axes.  This  orientation  can  also  be  described  by  the 
orientation  angles  (</>g,br  0?,b,  Yg,t)-  These  two  sets  of  angles  are  related  as  follows: 


tana  =  tan  6gb  cos  <j)gb  + 


tan  y/g'b  sin  4  , 


cos  0, 


g,b 


sin  p  =  cos  y/g  b  sin  0s  b  sin  <ps  b  -  sin  y/s  b  cos  $ 


'g,b 


g,b 


g,b 


rg* 


(5) 
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4.  Calculations  of  Grid  Positions 

Points  within  the  grid,  and  store  orientations  at  those  points,  are  defined  with  respect  to 
the  carriage  axis  system.  These  points  are  required  in  the  tunnel  axis  system  during 
testing,  and  can  be  converted  using  the  transforms  below. 

4.1  Transformations  between  Reference  Frames 

A  given  frame,  Fb,  is  related  to  another  frame,  Fv,  by  three  Euler  angles:  t//  (yaw),  9  (pitch) 
and  0  (roll).  To  transform  Fv  to  Fb,  these  rotations  are  applied  in  the  order  (y/,  9,  (f>),  and  to 
transform  Fb  to  Fv  the  rotations  are  applied  in  the  opposite  direction  and  in  reverse  order, 
(-0  -9,  -i//).  The  transformation  from  Fb  to  Fv  can  be  expressed  using  the  equation 

Fy  =  MvbFb  (6) 

where  the  transformation  matrix  Mv,b,  representing  the  rotations  (-$,  -6,  -y/),  is  [7]: 

Mv.b  =  Mi(->//).M2(-0).Mt(-0)=MTBJ,  (7) 

cos#cos^  sin 0 sin 0 cos y/ -  cos 0 sin y/  cos sin 6? cos ^  +  sm0smy/ 

=  cos<9sin^  sin  sin  6?  sin  ^  +  cos^cos^  cos  0  sin  0  sin  y/ -sin  0  cos  y/ 

-sin  0  sin  0  cos  0  cos  0  cos  0 

The  following  rotation  matrices  are  required  for  transformation  to  and  from  grid  points. 
Each  matrix  can  be  determined  by  substituting  the  appropriate  angles  into  Equation  7. 

Mc,b  =  Transformation  from  store  axes  to  carriage  axes 

-  order  of  rotations  (-0c,b,  -9c,b,  -  y/c,b ) 

Mac.c  =  Transformation  from  carriage  axes  to  aircraft  axes 

-  order  of  rotations  (-0ac,c,  -9ac,C/  -  yrac,c ) 

Mg,ac  =  Transformation  from  aircraft  axes  to  tunnel  axes 

-  order  of  rotations  {-0x,aC/  -9g:ltC/  -y/gAc) 

Transformation  matrices  can  be  combined  by  multiplication.  For  example,  converting 
frame  Fb  to  Fv  via  an  intermediate  frame  Fi  can  be  expressed  as 

Fv  =MvjM,bFb  (8) 

=  ^  V.B^B 

The  following  transformation  matrices  are  useful: 

Mg/C  =  Transformation  from  carriage  axes  to  tunnel  axes  =  Mg/acMaC/c  (9) 

Mg/b  =  Transformation  from  body  axes  to  tunnel  axes  ^ 

=  Mg,acMac,cMc,b  =  Mg,cMc,b 
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4.2  Transformation  of  Grid  Points  to  Tunnel  Axes 


The  grid  points  are  defined  in  terms  of  a  polar  angle  Op  and  radial  distance  r,  often 
normalised  by  the  store  diameter  D  to  produce  the  non-dimensional  radial  calibre,  r/D. 
The  grid  is  defined  in  carriage  axes  in  the  Yc-Zc  plane,  and  each  point  can  be  expressed  in 
Cartesian  coordinates  with  reference  to  carriage  axes  as 


=  0  (n) 

yc,h  =  r  sin  Op 
zc,b  =r  cosOP 

Using  the  transformation  from  carriage  axes  to  tunnel  axes,  the  grid  points  can  be 
expressed  in  tunnel  axes  relative  to  carriage  position  as 


(12) 


When  reducing  recorded  tunnel  data.  Equations  11  and  12  can  be  used  to  calculate  the  grid 
points  from  tunnel  axes,  such  that 


X 

(  \ 

Xc,b 

y g,b 

o 

U) 

II 

yc,b 

\zg’b  j 

gZc,bj 

(  \ 

r  \ 

Xc,b 

V 

yc,b 

=m;\ 

g,c 

KZc>bJ 

KZS’b) 

(13) 


f 

Op  =  arctan 

V 


(14) 


r  _  yj y?,b  +  Ze,6 

D~  D 


(15) 


Most  computer  languages  provide  a  function  atan2(y,x )  for  evaluating  the  inverse  tangent 
(arctan)  with  two  arguments,  as  shown  in  Equation  14.  This  function  produces  angles 
between  -n  and  n.  Most  languages  define  the  first  argument  as  the  numerator  and  the 
second  as  the  denominator.  Other  languages  including  Microsoft  Excel  have  the  order  of 
the  two  arguments  reversed,  so  care  should  be  taken  during  implementation. 


4.3  Transformation  of  Store  Orientations  to  Tunnel  Axes 


Store  orientation  angles  are  defined  with  respect  to  carriage  axes,  as  required  by  DSTOres 
and  STEME  software.  Sometimes  it  is  desirable  to  express  store  orientation  with  respect  to 
tunnel  axes  at  each  grid  point.  This  can  be  found  from  the  elements  of  the  transformation 
matrix  Mg,b  of  Equation  10,  where  Mg,b,  similar  to  Equation  7,  is  given  by 


UNCLASSIFIED 


11 


DSTO-TN-1339 


UNCLASSIFIED 


C0g 

,bcvg,b 

s^bseg,c¥g,-c^bs¥g. 

ctgj>s08j>cV'gjb  +  stgjbsV'i 

,bS^S,b 

gjb  +  C^gbCy/gt 

C<t>:i,bsegbs¥gl,  -  S(/)g  bCy/, 

- 

s^bcog. 

C0g,bC0g,b 

Mu 

m12 

Ma~ 

m21 

M22 

M23 

m3\ 

m32 

^33  _ 

g,b 

gJ> 


(16) 


For  brevity,  "C"  represents  the  trigonometric  function  cosine,  and  "S"  represents  sine  in 

Equation  16.  Mtj  represents  the  element  of  the  Mg,b  matrix  in  the  xh  row  and  f1  column. 
From  the  elements  of  Mg/b  it  is  easy  to  obtain 


6 '  b  =arcsin(-M31) 


=  arctan 


"XT 

VM33  j 


f  m21a 

¥  b  =  arctan 


Mn 


(17) 

(18) 

(19) 


These  equations  are  valid  for  all  0,;h  and  ¥grb  within  the  range  specified  in  Equation  4. 
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5.  Determination  of  Aerodynamic  Loads 

Aerodynamic  loads  at  particular  store  positions  and  orientations  are  determined  by  post¬ 
processing  data  gathered  from  wind  tunnel  tests.  This  data  is  required  for  store  trajectory 
prediction  using  programs  such  as  DSTOres  and  STEME.  Section  5.1  describes  the  process 
by  which  raw  strain  gauge  balance  measurements  can  be  transformed  into  forces  and  vice 
versa.  Sections  5.2  to  5.13  detail  the  full  post-processing  procedure  used  to  produce  the 
required  data  for  store  trajectory  prediction.  A  diagram  of  the  post-processing  procedure 
is  shown  in  Figure  7. 


Figure  7  Post-processing  procedure  flow  chart 
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5.1  Reduction  of  Strain  Gauge  Balance  Outputs 


Each  of  the  output  channels  of  a  six-component  strain  gauge  balance  is  represented  by  a 
27-term  second  order  polynomial  equation  of  the  form  [8,  §4.2]: 


R,  =  Ct,H,  +  Ci2H2  +  Ci3H3  +  Ci4H4  +  CiiHi  +  Ci6H6 

+  CmH f  +  Ci22H2  +  Cu33H]  +  CiMH\  +  Ci55H25  +  CmH\ 

+  Cil2HlH2  +  CiX3HxH3  +  CiUHlH4  +  C^H,  +  Cil6HxH6 
+  Ci23H2H3  +  Ci24H2H4  +  Ci25H2H5  +  Ci26H2H6 
+  Ci34H3H4  +  Ci35H3H5  +  C.36H3H6 
+  CiA5H4H5+CiA6H4H6 

+  c,56h5h6 


i  =  1 . 6 


(20) 


where  the  C's  are  the  calibration  coefficients  and  Hi  are  component  loads. 

The  set  of  six  equations  represented  by  Equation  20  can  be  expressed  in  matrix  form: 

R  =  [C1]H  +  [C2]H*  (21) 

where 

R  is  the  vector  of  strain  gauge  balance  outputs 

[Cl]  is  the  matrix  of  linear  calibration  coefficients 

[C2]  is  the  matrix  of  non-linear  calibration  coefficients 
H  is  the  vector  of  component  loads 

H*  is  the  vector  of  squared  and  cross  products  of  component  loads 

The  calibration  matrices  are  normalised  by  pre-multiplying  each  by  a  diagonal  matrix  [D] 
composed  of  the  reciprocal  of  the  diagonal  elements  of  [Cl]  giving: 

[XI]  =  [D]  [Cl]  (22) 

[X2\  =  [D]  [C2] 

Equation  21  can  therefore  be  written  as: 

[D]R  =  [D][C1]  H  +  [D][C2]  H  *  (23) 

or 

H'=[X1]H  +  [X2]H*  (24) 

where  H'  is  a  vector  of  approximate  component  loads.  Solving  for  H  gives: 

H  =  [Xl]~x  H'~  [Xl]~l  [X2]  H  *  (25) 

An  iterative  scheme  is  used  to  determine  H  according  to  Equation  25. 

A  first  approximation  to  H  is  obtained  by  ignoring  the  non-linear  terms: 
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H,  =  [Xiy1  H'  (26) 

This  is  then  used  to  obtain  a  first  approximation  to  the  non-linear  term,  Ai,  thus: 

M  =  [XI]- 1  [X2]  Hi*  (27) 

where  the  elements  of  the  vector  Hi*  are  composed  of  the  squares  and  cross- 
product  of  the  elements  of  the  approximate  vector  Hi. 

A  second  approximation  H2  is  then  obtained  from: 

H2  =  H,-Ai  (28) 

The  nth  approximation  H„  is  obtained  iteratively  by  means  of  the  generalised  equations: 

An-1  =  [Xl]-1  [X2\  Hn-1*  (29) 

Hn  —  Hj  An-1 

For  most  balances  used  in  DSTO,  the  values  of  H  typically  converge  to  within  ±0.0001  N  or 
Nm  after  2  iterations. 

5.2  Zero  Corrections  for  Strain  Gauge  Balance  Outputs 

Zero  readings  are  taken  at  the  beginning  and  end  of  each  test  run.  For  each  zero  reading, 
the  model  orientation  is  set  to  zero  with  respect  to  tunnel  axes  (y/g,b  =  9g,b  =  <j)g,b  =  0°).  These 
initial  and  final  readings  are  averaged  to  produce  the  wind-off  zeros,  R0.  (  i  =  1,  •  •  •  ,6 ). 

5.3  Calculation  of  Buoyant  Tares 

The  tare  loads  applied  to  the  strain  gauge  balance  by  the  metric  mass  of  the  store  model  at 
9g,b  =  0°  and  <f>g,b  =  0°  are  computed  according  to  the  following  equations  [8,  §4.3]: 

-  kx  sin  6gb  (30) 

kY  cos  0gb  sin  (j)g  b 
kz  cos  9gb  cos  <f>gb 

k-Mxi  cos  0gb  cos  (f>gb  —  kMXX  cos  6g  b  sin  (f>g  b 
-  kMY2  sin  0g  b  -  kMY[  cos  6g  b  cos  </>g  b 
k\tz2  cos  6g  b  sin  (f>g  h  +  kMZX  sin  6g  b 

where  kx,  ky,  kz,  fevrxi,  /cmx2,  /cmyi,  kMYi,  kMzi,  fcuz2  are  the  tare  weight  coefficients 
determined  for  each  model  configuration  prior  to  testing. 

The  readings  of  the  six  strain  gauge  balance  outputs,  Rr.  ( i  =  1,  •  •  •  ,6 ),  corresponding  to 
these  loads  are  calculated  from  the  strain  gauge  balance  calibration  matrices  according  to 
Equation  21. 

A  set  of  buoyant  zero  readings,  R0.  ( i  =  1, •  •  •  ,6 ),  are  obtained  by  subtracting  these  readings 
from  the  wind-off  zeros,  R'0j  ( /  =  1, •  •  •  ,6 ),  thus: 
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*0/  =  *0/  —  -R-Ti  ?  ^  =  1?  *  *  *  (31) 

These  buoyant  zeros  are  used  to  correct  wind-on  strain  gauge  balance  outputs. 

5.4  Reduction  of  Strain  Gauge  Balance  Outputs  to  Gross  Loads 

Wind-on  measurements  produce  six  raw  component  measurements  from  the  strain  gauge 
balance,  (R. ).  The  buoyant  zeros  are  subtracted  from  these  raw  measurements  to  give 

absolute  strain  gauge  balance  readings,  (R.  ): 

Ri  =Rl-R0i  >  i  =  (32) 

These  absolute  readings  are  converted  into  absolute  loads  according  to  the  process 
described  in  Section  5.1.  These  absolute  loads  are  termed  the  strain  gauge  balance  gross 
loads  (Hg). 

5.5  Calculation  of  Sting-Balance  Deflection 

Under  load,  the  sting  and  strain  gauge  balance  structure  deflects.  This  deflection  can 
generally  be  considered  in  combination.  The  deflections  about  the  balance  moment  centre 
can  be  determined  using  the  following  methodology. 

The  components  of  the  sting-balance  deflection  in  balance  axes  are  the  displacements  x  ,  y', 
and  z'  and  the  angular  deflections  v,  77  and  %.  These  are  related  to  the  applied  loads,  Hg ,  by 
the  stiffness  matrix  K,  which  consists  of  10  experimentally-derived  coefficients,  ki,  kz,  . . ., 
kw : 

^10  0  0  0  0  0  Fxbalg 

0  0  0  0  k9  Fybalg 

0  0  k6  0  k,  0  Fzbalg 

0  0  0  ^00  Mxbalg 

0  0  k2  0  k3  0  Mybalg 

_  0  k4  0  0  0  k5\[MZBALG 

5.6  Transformation  of  Store  Position  and  Orientation 

The  model  position  and  orientation  relies  on  the  readings  from  the  store  support  (SS) 
system.  Corrections  need  to  be  applied  due  to  the  sting-balance  deflections  (see  Section 
5.5)  and  any  misalignment  between  the  balance  and  model  body  axes  (see  Section  5.8).  The 
store  support  system  computes  and  outputs  parameters  in  tunnel  gravity  axes  based  on 
feedback  from  angular  sensors  and  linear  displacement  transducers  (LDTs).  The  axial  (xgrSS) 
and  vertical  (zg/SS)  positions  are  measured  by  dedicated  LDT's.  The  pitch  angle  pitch  (#g/SS) 
is  measured  by  a  dedicated  encoder.  Roll  angle  (^/SS)  is  measured  by  a  resolver.  The  lateral 
position  (y?/S5)  and  yaw  angle  ( i//g/SS)  are  calculated  from  the  measurements  of  a  pair  of  yaw 
encoders. 
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For  small  misalignment  angles  between  the  strain  gauge  balance  and  the  model,  the  store 
support  yaw,  pitch,  and  roll  offsets  are  nulled  during  initial  model  setup.  Strain  gauge 
balance  misalignment  is  subsequently  disregarded  in  the  calculation  of  the  model  position 
and  orientation.  In  cases  where  there  is  a  large  misalignment  between  the  strain  gauge 
balance  and  either  the  model  or  store  support  system,  care  must  be  taken  to  properly 
account  for  the  misalignment  during  setup  and  post-processing.  These  cases  are  beyond 
the  scope  of  this  document. 


The  strain  gauge  balance  deflection  angles  are  not  Euler  angles,  but  as  they  are  small  they 
can  be  approximated  by  Euler  roll,  pitch  and  yaw  angles  [9,  p.108],  resulting  in  the  Euler 
rotational  transformation  matrix  MSS/b,  (v,  77  and  %).  The  Euler  rotational  transformation 
matrix  from  the  store  support  system  orientation  to  tunnel  gravity  axes  is  Mg/SS/  (fass,  0g/SS 
and  i//g,ss)-  The  model  position  in  tunnel  gravity  axes  is  given  by: 


Xg,ss 

~xr 

y8 

= 

y  g,ss 

+  Mg,s 

y 

_v 

_Zg,ss  _ 

z' 

(34) 


The  model  orientation  with  respect  to  tunnel  gravity  axes  can  be  derived,  as  described  in 
Section  4.3,  from  the  elements  of  the  Euler  rotational  transformation  matrix  Mg:b,  which  is 
calculated  as: 


Mg,b  =Mg>ssMss,b 


(35) 


The  location  of  the  store  model  in  tunnel  axes  relative  to  carriage  position,  (xg,t,  yg,b,  zg,b),  can 
be  found  by  subtracting  the  carriage  position  in  tunnel  axes,  (xg/Cr  yg,c,  zg,c),  from  the  model 
position  in  tunnel  axes,  (xg,  yg,  zg): 


x  ,-x  -x 

g,b  g  g,c 


(36) 


yg,b =yg~y 


g,C 


Z  u  —  z  ~  Z 

g,b  g  g,c 


The  store  grid  coordinates  r/D  and  Op  can  be  calculated  using  Equations  13, 14  and  15. 


The  aircraft  angle-of-attack  ( 6g,ac )  is  derived  directly  from  the  parent  aircraft  turntable 
encoder  reading,  and  no  data  reduction  is  required  from  the  raw  tunnel  output  file.  The 
carriage  offsets  ( 0ac,c  and  ^flc,c)  are  known  in  advance  from  model  geometry  and  are  not 
measured  during  testing.  These  angles  can  be  used  to  form  the  rotation  matrix  Mg/C 
(Equations  7,  9).  The  rotation  matrix  Mg/b  is  calculated  in  Equation  35.  Equation  10  can  then 
be  rearranged  to  get  Lc,b, 


M  h  =  M  M 

c,b  g,c  g,s 


(37) 


The  elements  of  Mc,b  can  be  used  to  calculate  the  store  orientation  angles  in  carriage  axes 
{</>c,br  0C/b  and  y/C/b),  in  the  same  way  that  Equations  17, 18  and  19  are  used  to  calculate  store 
orientation  in  tunnel  axes  from  Mg,b. 
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5.7  Calculation  of  Aerodynamic  Loads  in  Balance  Axes 

The  tare  loads  at  the  model  attitude  associated  with  the  wind-on  measurement  (Hr)  are 
computed  using  Equation  30.  Note  that  model  attitude  must  be  relative  to  tunnel  axes,  as 
calculated  from  the  elements  of  the  matrix  Mg/b  in  Equation  35  (i.e.  model  orientation 
corrected  for  deflections). 


The  aerodynamic  loads  in  balance  axes  (H)  are  then  calculated  by  subtracting  tare  loads 
from  gross  loads: 


5.8  Transformation  of  Aerodynamic  Loads  from  Balance  to  Body  Axes 

To  account  for  any  small  rotational  misalignment  between  the  balance  axes  and  the  body 
axes  of  the  model,  the  transformation  matrix  used  can  be  expressed  in  terms  of  the 
following  offset  angles: 

gi  —  Angle  between  the  projection  of  the  model  Xb  axis  in  the  Xbai-Ybai  plane  and 
the  Xbai  axis,  positive  nose  to  starboard. 

£2  —  Angle  between  the  model  Xb  axis  and  its  projection  on  the  Xbai-Ybai  plane, 
positive  nose  up. 

<^3  —  Angle  between  the  model  Zb  axis  and  its  projection  on  the  Xbai-Zbai  plane, 
positive  clockwise  looking  forward. 


Note  that  each  of  these  angles  rotate  about  the  undisturbed  axis,  as  distinct  from  a  set  of 
Euler  angles.  The  rotational  transformation  matrix  from  balance  axes  to  store  body  axes, 
derived  in  Appendix  B.l,  is  given  by: 


M 


b,bal 


cos  cos  <^2  sin  ^  cos 

sin  sin  ^3  -  A  cos  ^  sin  ^  cos2  i;2  A 

_______________________________ 

B 

cos  sin  £  cos2  £  sin  £  +  A  sin  &  ^  e 

sm  u 

B  3 


-  sin 

cos  <^2  (cos  £  sin  <^3  +  A  sin  ^  sin  ^2 ) 
B 

cos  £2  (-  sin  ^  sin  %2  sin  ^3  +  A  cos  ^ ) 
B 


(39) 

where 
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A  =  ^/cos2  -  sin2  ^  cos2  <^2 

and 

B  =  cos2  ^  cos2  <^2  +  sin2  £2 


The  transformation  of  the  aerodynamic  forces  from  balance  axes  to  body  axes  is  therefore 
given  by: 

F=MwFbal  (40) 


where 


~FX~ 

F 

1  XBAL 

F  = 

Fy 

and 

F  = 

1  BAL 

F 

1  YBAL 

Fz_ 

F 

|A  ZBAL  J 

are  the  vectors  of  aerodynamic  force  in  body  and  balance  axes  respectively. 

To  account  for  the  translational  misalignment,  let  (pcb,  yb,  Zb)  be  the  location  of  the  model 
moment  reference  centre  (MRC,  origin  of  the  body  axes  system)  with  respect  to  the  balance 
moment  centre  (BMC,  origin  of  the  balance  axes  system).  Then  the  transformation  of  the 
aerodynamic  moments  from  balance  axes  to  body  axes  are  given  by: 


"  0 

zb 

-yb 

M  ~  Mb,balMBAL  + 

0 

xb 

F 

_  yb 

~Xb 

0 

where 


~MX~ 

M  XBAL 

M  = 

My 

,  and 

M BAL  ~ 

^  YBAL 

_M  z  _ 

ZBAL  _ 

(41) 


are  the  vectors  of  aerodynamic  moment  in  body  and  balance  axes  respectively. 

Since  the  rotations  are  applied  first,  followed  by  translation,  Xb,  yb,  Zb  must  be  measured  in 
body  axes. 


5.9  Transformation  to  Body  Axes  with  Origin  at  Centre  of  Gravity 

For  use  with  legacy  versions  of  DSTOres,  forces  and  moments  must  be  defined  about  the 
centre  of  gravity  (c.g.)  of  the  real  store.  This  will  be  orthogonal  to  the  model  axes  system 
but  may  have  its  origin  at  a  different  position  than  the  nominal  MRC.  To  account  for  this 
translational  misalignment,  let  (xcg/  ycg/  zcg )  be  the  location  of  the  model  c.g.  with  respect  to 
the  MRC.  The  forces  applied  at  the  c.g.  do  not  change,  and  the  translation  of  the  model 
aerodynamic  moments  about  the  c.g.  are  given  by: 
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Mc.g=M  + 


0 

•U* 


eg 

0 


■y, 


eg 


"eg 

0 


(42) 


5.10  Calculation  of  Flow  Conditions 


The  test  section  flow  conditions  are  calculated  from  the  reference  transducers,  settling 
chamber  total  pressure  (Pf),  total  temperature  (Tf),  and  plenum  static  pressure  (Ps),  using 
isentropic  flow  equations  as  detailed  below. 


5.10.1  Mach  number  (M) 


M  = 


2 


P^ 

~Pt 


7 


where  y=1.4  is  the  ratio  of  specific  heats  for  air 


5.10.2  Dynamic  pressure  (q) 


(43) 


q  =  —PM2  =  0.1PM2 
2 

5.10.3  Reynolds  number  (Re) 


P 

7  PT 


(1+2M) 

Md  1 1.4  f, 

“  q«r,  (1  +  0.2M2/5 

where  p  is  the  flow  density 
V  is  the  flow  velocity 
d  is  the  model  reference  dimension 


(44) 


(45) 
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R  is  the  gas  constant 

Ts  is  the  static  temperature  (calculated  from  the  total  temperature 
and  Mach  number) 

/j  is  the  kinematic  viscosity  (calculated  from  Sutherland's  law) 


5.11  Conversion  to  Coefficient  Form 

The  aerodynamic  force  coefficients,  (Cx,  Cy,  Cz),  and  corresponding  moment  coefficients, 
(C/,  Cm,  Cn ),  acting  on  the  model  can  be  calculated  from  the  aerodynamic  forces  and 
moments  in  body  axes: 


F 

r  —  x  • 

V- '  y  —  / 

qS 

r  =^JL- 

'  qSd 


c  -A. 

—  / 


Q=- 


c  c  =h4^- 

m  qSd  "  qSd 


where  q  is  the  dynamic  pressure,  S  is  the  model  reference  area,  and  d  is  the  model 
reference  dimension. 

The  dynamic  pressure  is  calculated  from  the  reference  pressure  transducers  in  accordance 
with  Equation  44. 

5.12  Transformation  to  Wind  Axes 

Aerodynamic  force  and  moment  coefficients  can  be  converted  from  store  body  axes  to 
wind  axes  using  the  transformation  matrix: 

cos  a  cos  (3  sin  (3  sin  a  cos  (3  (47) 

M wb=  -cos  a  sin  (3  cos  (3  -sin  a  sin  /3 
-sin  a  0  cos  a 

The  aerodynamic  angles  a  and  /?  are  calculated  using  Equation  5. 


5.13  Transformation  to  Missile  Axes 


Aerodynamic  force  and  moment  coefficients  can  be  converted  from  store  body  axes  to 
missile  axes  (non-rolling  axes)  using  the  transformation  matrix: 

"1  0  0  1  (48) 

MPjb=  0  cos^c5  -  sin  (f>cs 
0  sin  (j)cs  cos  4,  _ 
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6.  Conclusion 


This  document  has  presented  the  algorithms  required  to  prepare  test  plans  for  transonic 
wind  tunnel  grid  testing  and  to  reduce  wind  tunnel  grid  data  for  store  release  analysis. 
Elements  of  the  procedures  presented  may  also  be  amenable  to  other  wind  tunnel  tests, 
and  the  algorithms  have  been  presented  in  general  form  to  assist  with  this.  The 
implementation  of  these  algorithms,  while  not  covered  in  this  document,  is  an  important 
task  that  should  be  thoroughly  verified  and  validated  to  prevent  avoidable  errors  and 
uncertainty  from  affecting  the  final  results.  The  example  calculations  in  Appendix  C  can  be 
used  to  validate  software  implementations  of  these  algorithms. 
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Appendix  A:  The  DSTO  Transonic  Wind  Tunnel 


DSTO  maintains  and  operates  several  experimental  aerodynamic  test  facilities  including  a 
transonic  wind  tunnel  (TWT).  The  TWT  was  commissioned  by  the  US  supplier.  Aero 
Systems  Engineering  (ASE).  It  is  the  only  facility  of  its  kind  in  Australia,  and  enables 
aerodynamic  research  to  be  carried  out  for  the  Australian  Defence  Force  (ADF).  The 
facility  operates  in  the  transonic  speed  range  where  both  viscosity  and  compressibility 
effects  are  important  and  are  difficult  to  model  using  other  means.  This  facility  is  a  critical 
element  of  Australia's  indigenous  aircraft-stores  combatibility  clearance  capability. 


The  TWT  is  a  closed  return-circuit  continuous  flow  tunnel.  A  schematic  drawing  of  the 
aerodynamic  circuit  and  main  components  of  the  facility  is  shown  in  Figure  Al.  The 
facility  has  a  nominal  Mach  number  range  from  0.3  to  1.2  continuously  variable,  and  Mach 
1.4  with  a  fixed  nozzle.  The  total  pressure  can  be  varied  from  30  kPa  to  200  kPa.  Flow  is 
generated  by  a  two-stage  axial  flow  compressor  powered  by  a  5.3  MW  variable  speed 
electric  motor.  Each  stage  comprises  of  30  rotor  blades  and  the  fan  speed  can  be  driven  to 
1800  rpm.  The  blade  positions  are  manually  adjustable  however  they  were  fixed  at 
predefined  positions  during  commissioning. 


Figure  Al  The  DSTO  transonic  wind  tunnel  aerodynamic  circuit 

The  tunnel  is  equipped  with  a  Plenum  Evacuation  System  (PES)  to  pressurise  and 
depressurise  the  tunnel  circuit,  and  to  improve  flow  quality.  The  PES  consists  of  a  17  blade 
centrifugal  compressor  driven  by  a  2.6  megawatt  motor  and  operates  at  a  constant  speed 
of  10840  rpm.  This  extracts  air  from  the  plenum,  compresses  it  by  a  factor  of 
approximately  2.9  and  then  injects  the  compressed  air  downstream  of  the  first  corner.  The 
effect  that  this  has  on  the  test  section  flow  is  equivalent  to  introducing  a  steady  increase  in 
cross-sectional  area  in  the  downstream  direction. 
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A  water-cooled  heat  exchanger  in  the  settling  chamber  is  used  to  control  the  air 
temperature  in  the  test  section.  Turbulence  conditioning  screens  and  a  16:1  contraction 
provide  good  quality  flow. 

Three  independent  model  support  systems  are  provided: 

•  main  model  support  for  freestream  tests  of  a  full  model,  providing  pitch  and  roll 
capability 

•  store  model  support  to  position  a  store  model  below  a  parent  model  for  store 
release  tests,  with  full  six  degree  of  freedom  movement 

•  sidewall  turntable  model  support  for  half-model  tests,  providing  pitch  rotation 

The  test  section  is  nominally  2700  mm  long  and  has  a  cross  section  of  806  mm  square.  The 
test  section-nozzle  interface  is  at  station  3505.  The  test  section  model  support  interface  is  at 
station  6205.  The  high  flow  quality  region  extends  ±400  mm  from  the  nominal  model 
centre  at  station  5257.  The  test  section  is  shown  in  Figure  A2  together  with  the  other 
components  which  form  the  test  leg. 

The  test  section  sidewalls  are  interchangeable  and  can  be  either  solid  or  slotted.  Both  the 
top  and  bottom  wall  angles  were  fixed  to  +20  minutes  each  (divergent)  during 
commissioning.  However  they  are  manually  adjustable  to  ±  30  minutes  each  of  arc.  Each 
slotted  wall  consists  of  six  slots  per  wall  with  a  porosity  of  4.97%. 


UNCLASSIFIED 


25 


UNCLASSIFIED 

DSTO-TN-1339 


UNCLASSIFIED 

DSTO-TN-1339 

Appendix  B:  Transformation  Derivations 

B.l.  Derivation  of  the  Rotational  Transformation  from  Balance  Axes  to 
Body  Axes 

The  following  offset  angles  describe  the  rotational  misalignment  between  the  balance  axes 
and  the  body  axes: 

£1  —  Angle  between  the  projection  of  the  model  Xb  axis  in  the  Xbai-Ybai  plane  and 
the  Xbai  axis,  positive  nose  to  starboard 

—  Angle  between  the  model  Xb  axis  and  its  projection  on  the  Xbai-Ybai  plane, 
positive  nose  up 

<^3  —  Angle  between  the  model  Zb  axis  and  its  projection  on  the  Xbai-Zbai  plane, 
positive  clockwise  looking  forward 

These  angles  are  not  Euler  angles,  but  their  transformation  matrix  can  be  derived  by 
considering  three  Euler  angle  rotations: 

Mi  —  Rotate  by  §  around  the  model  Zb  axis. 

M2  —  Rotate  by  &  around  the  model  Yb  axis. 

M3  —  Rotate  by  (j)  around  the  model  Xb  axis. 

F  =  M3(m2^2W^l)FBAL  =  MbMlFBAL 


Elements  An,  An  and  An  are  in  terms  of  £1  and  ^2  only,  and  are  therefore  derived.  For  the 
rest  of  the  elements,  a  relationship  must  now  be  found  between  </)  and  This  can  be  found 
by  considering  the  components  of  Zb  in  balance  axes,  illustrated  in  Figure  Bl. 
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Figure  B1  Components  of  Zb  in  Balance  Axes 
From  the  Ybai  component,  we  have 

sin(-£)  =  ^32 

-  sin  £3  =  -  cos  <fjj  sin  (/>  +  sin  sin  cos  (/> 

From  the  projection  of  Zb  on  the  Xbai-Zbai  plane,  we  have 
cos(-^3)=a/4T4 

cos2  =  (sin  ^  sin^  +  cos^  sin^  cos^)2  +(cos^2  cos^)2 
On  expansion  and  rearrangement,  it  can  be  shown  that 

cos2  <J3  -  sin2  £  cos2  <^2  =  (cos ^  cos^  +  sin^  sin  £2  sin^)2 


For  convenience,  let 

A  =  -^/cos2  <^3  sin2  ^  cos2  £2  =  cos  ^  cos  (j>  +  sin  ^  sin  <^2  sin 


Hence, 

A22  =  cos^  cos^  +  sin^  sin i;2  sin^  =  ^/cos2  <^3  -  sin2  ^  cos2  £2  =  ^4 

B.1.1  Deriving  A22 

Using  Equations  B2  and  B3,  A21  can  be  expressed  in  terms  of  £1,  ^2  and  ^3. 
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A2l  =  -  sin  cos  (j)  +  cos  sin  sin  (j) 

_  (-  sin  ^  cos  (j)  +  cos  ^  sin  sin  ^)(cos2  ^  cos2  +  sin2  ^2 ) 
cos2  cos2  ^2  +  sin2 


Let 


S  =  cos2  cos2  +  sin2 

i21xg  =  -sin^j  cos2  cos2  cos^-sin<^  sin2  £2  cos^ 

+  cos  sin  <^2  cos2  ^  cos2  ^2  sin  (f>  +  cos  ^  sin  ^2  sin2  E,2  sin  </> 

=  sin  E,2  (cos  sin  (j)  -  sin  ^  sin  ^2  cos  (f)-  cos  ^  sin  E,2  cos2  <^2  sin  <j> 
-  sin  cos2  ^  cos2  ^2  cos  ^  + cos  ^  sin^2  cos2  ^  cos2  ^2  sin^ 


Using  Equation  B2 


A2lxB  =  sin  E,2  sin  (^3  -cos^  sin<^2  cos2  sin  sin  ^  cos2  ^  cos2  cos^ 
+  cos^j  sin  cos2  ^  cos2  sin^J 

=  sin  t;2  sin  £3  -  cos  sin  ^  cos2  (sin  ^  sin  E,2  sin  (j)  +  cos  ^  cos  <f) 
Using  Equation  B3 


An  — 


sin  sin  £3  -  A  cos  sin  cos2 
B 


B.1.2  Deriving  A32 

Using  Equations  B2,  B3  and  B4,  element  A31  can  be  expressed  in  terms  of  £1,  6  and  £3. 
A31  =  sin  sin  (f> + cos  sin  cos  (f) 

_  (sin^j  sin  (/)  +  cos  ^  sin  £2  cos^)(cos2  ^  cos2  +sin2  £2) 
cos2  ^  cos2  £2  +  sin2 

Using  Equation  B4 

A3lxB  -  sin^  cos2  £  cos2  ^2  sin  ^  + sin  sin2  sin^ 

+  cos^j  sin  cos2  ^  cos2  cos  ^  + cos  ^  sin^2  sin2  cos^ 

=  cos  sin  cos2  (cos  ^  sin  (/)  -  sin  ^  sin  <^2  cos  (j)) 

+  sin  (sin  ^  sin  sin  (j)  +  cos  ^  cos2  cos  (j)  +  cos  ^  sin2  cos 

Using  Equation  B2 

A3lxB  =  cos  sin  ^  cos2  sin 

+  sin  (sin  ^  sin  sin  (j)  +  cos  ^  cos2  cos  (j)  +  cos  sin2  cos 

=  cos^j  sin^j  cos2  sin^3  +  sin£2(sin£1  sin£,  sin (f>  +  cos ^  cos^) 


(B4) 
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Using  Equation  B3 


'31 


cos  sin  ^  cos2  sin  +  A  sin  A2 
5 


B.1.3  Deriving  A23 

Using  Equations  B2,  B3  and  B4,  element  A23  can  be  expressed  in  terms  of  ca,  A  and  A- 

A23  =  cos  A2  sin  (j) 

_  (cos  A2  sin^)(cos2  A  cos2  A2  +  sin2  £2) 
cos2  A  cos2  A2  +  sin2 


Using  Equation  B4 


A23xB  =  cos  ^2  sin  ^(cos2  A  (l  -  sin2  A2 )+  sin2  A2 ) 

cos  A  (cos  A  sin  (j)  -  sin  A  sin  £2  cos  (f)  ' 

v+  sin  A  sin  A2  (cos  A  cos  (j)  +  sin  A  sin  A2  sin  ^)y 


=  cos  A 


Using  Equations  B2  and  B3 


A 


cos  A2  (cos  sin  ^  ^  sin  A  sin  A2 ) 


■23 


5 


B.1.4  Deriving  A33 

Using  Equations  B2,  B3  and  B4,  element  A33  can  be  expressed  in  terms  of  A/  A  and  A- 

A33  =  cosA2  cos^ 

_  (cos  cos^)(cos2  A  cos2  A2  +  sin2  £2) 
cos2  Ai  cos2  ^2  +  sin2 


Using  Equation  B4 


A33xB  =  cosA2  cos^(cos2  ^(l-sin2  ^2)+sin2  A2) 
^cos  A  (cos  A  cos  </>  +  sin  A  sin  £2  sin 


=  cosA2 


-  sin  A  sin  g2  (cos  A  sin  (f)  -  sin  A  sin  g2  cos 


Using  Equations  B2  and  B3 


^33  _ 


cos  ^2  (-  sin  ^  sin  A2  sin  +  A  cos  £ ) 


5 
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Appendix  C:  Sample  Calculations 

The  following  examples  are  presented  to  show  the  processes  required  to  prepare  for  grid 
testing  and  to  reduce  data  from  grid  testing  in  the  DSTO  transonic  wind  tunnel.  These 
examples  can  also  be  used  to  validate  software  developed  for  this  purpose.  Note  that  the 
numbers  here  are  not  necessarily  realistic,  and  precision  is  retained  for  validation 
purposes,  not  as  an  indication  of  uncertainty. 

C.l.  Converting  a  Test  Point  to  Tunnel  Axes 

The  test  point  is  defined  inTable  Cl. 


Table  Cl  Test  point  definition 


Variable 

Symbol 

Value 

Units 

Store  diameter 

D 

0.02457 

m 

Grid  radial  calibre  position 

r/D 

3.5 

ND 

Grid  lateral  polar  angle 

0p 

-20 

degrees 

Parent  aircraft  angle-of-attack 

9%,ac 

6 

degrees 

Carriage  axes  pitch  offset 

Oac,c 

-3 

degrees 

Carriage  axes  roll  offset 

(frac,c 

5.5 

degrees 

Store  yaw  angle  (carriage  axes) 

Yc,b 

5 

degrees 

Store  pitch  angle  (carriage  axes) 

9c,b 

-15 

degrees 

Store  roll  angle  (carriage  axes) 

(j)c,b 

10 

degrees 

C.1.1 


Rotation  Matrices 


The  rotation  matrices  can  be  calculated  using  Equation  7.  Note  that  all  angles  not  defined 
in  Table  Cl  are  0.  The  rotation  matrices  are: 


M 


g,ac 


{-0g,J 

cos  egac  o 

Sin0g,ac 

0  1 

0 

_~Sin0g,ac  0 

COS0g,ac_ 

'  0.9945219 

0  0.1045285] 

0 

1  0 

-0.1045285 

0  0.9945219] 
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C°S  0ac,c 

sin  fiac.c  sin  9acx 

cos (jfacc  sin  0acc 

0 

cos  tc,c 

-sin  tc,c 

-sin  dac,c 

sin  fiac.c  cos^C;C 

cos  <t>aC'C  cos0acc 

0.9986295 

-0.0050162 

-0.0520950" 

0 

0.9953962 

-0.0958458 

0.0523360 

0.0957144 

0.9940320 

M,t  =M2(-Vcl\M2  (-  6,„  ]M,  (-  ti ) 

cos  9cb  cos  Xf/C  b  sin  (/>c  b  sin  0ch  cos  y/cb  -  cos  </>c  b  sin  \f/c  b 
=  cos  0c  b  sin  y/c  b  sin  (/>c  b  sin  0c  b  sin  y/c  b  +  cos  </>c  b  cos  y/c  b 

-sin  0c,b  sin^c6  cos6cb 

cos  (f)c  b  sin  0c  h  cos  y/cjb  +  sin  <f>c  b  sin  y/cjb 
cos  (j)c  b  sin  0c  h  sin^c  6  -sin^e6  cos^c  6 
cos^fecos6»c6 

"0.9622502  -0.1306041  -0.2387826" 

=  0.0841860  0.9771432  -0.1952023 

0.2588190  0.1677313  0.9512513 

0.0520950 " 
-0.0958458 
0.9940320 


=  Mg,acMuc,c 


0.9986295  0.0050162 

0  0.9953962 

-0.0523360  0.0957144 


Mg,b  =  Mg^cMaccMct)  =  M  gcMcb 


0.9748369  -0.1167856  -0.1898791 
0.0589917  0.9565683  -0.2854770 

0.2149719  0.2670922  0.9393874 
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C.1.2  Grid  Position  in  Tunnel  Axes 
The  grid  position  in  carriage  axes  is: 

xc,b=  0m 

ycb=rsin&P 

= -0.0294120  m 
zcJ>=r  cos&P 
=  0.0808089  m 


where  r  =  r/D*D  =  0.085995  m. 


The  grid  position  in  tunnel  axes  is  calculated  as  follows: 


f  \ 

/  \ 

Xcjb 

y  g,b 

0 

bo 

11 

yc,b 

\zg’b ) 

KZ°’b  J 

'  0.0040622  ' 
-0.0370218 
v  0.0775114  J 


m 


C.1.3  Store  Orientation  in  Tunnel  Axes 


The  store  orientation  in  tunnel  axes  is  found  from  the  elements  of  Mg/b  in  accordance  with 
Equations  17, 18  and  19: 

MgJ>  (3,  l)=sin0gb  =0.2149719 


0gb  =-12.4139° 


Mg  b  (3,2)  =  sin  </>gb  cos  6g  b  =  0.2670922 
Mgb( 3,3)  =  cos 4,  cos Ggb  =  0.9393874 


^g,b 


( 

=  arctan 

v 


=  15.8718° 


M  gh{\,\)  =  co$0gh  cos  y/g>b  =  0.9748369 
Mgb{2,\)  =  cos0gb  sin  y/gb  =  0.0589917 


f 

y/gb  =  arctan 

v 


Mjy if 


=  3.4630° 
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C.2.  Data  Reduction  for  a  Test  Point  Measurement 

The  general  test  data  for  this  example  are  listed  in  Table  C2,  strain  gauge  balance 
coefficients  are  listed  in  Table  C3,  and  data  specific  to  this  particular  test  point  are  listed  in 
Table  C4.  Strain  gauge  balance  calibration  data  is  given  in  Equations  Cl  and  C2. 


Table  C2  Example  test  data 


Variable 

Symbol 

Value 

Units 

Store  diameter 

D 

0.02457 

m 

Store  reference  area 

Sref 

0.000474 

m2 

Location  of  MRC  w.r.t  BMC 

Smrc,bmc 

[-0.0025,  0,  0] 

m 

Zero  pitch  reading,  averaged 

00 

0.00206301998705 

degrees 

Zero  roll  reading,  averaged 

(j)0 

0.010111324 

degrees 

Gross  strain  gauge  balance  zero 
readings,  averaged 

R'o 

[0, 

0.000180001294, 

0.000769806525, 

0.000359504257, 

-0.000335901155, 

0.000030521535] 

mV/V 

Table  C3  Strain  gauge  balance  coefficients 


Coefficient 

Value 

Units 

Deflection  coefficients 

fa 

1.4689 

degrees/ Nm 

k2 

-0.0144 

degrees/ N 

fa 

0.4378 

degrees/ Nm 

h 

0.01447 

degrees/ N 

ks 

0.4515 

degrees/ Nm 

h 

0.03295 

mm/N 

k7 

-0.2403 

mm/Nm 

Tare  coefficients 

kx 

0 

N 

kY 

3.7399786 

N 

kz 

3.7844698 

N 

k]VLXl 

0.0084307101 

Nm 

k]VLX2 

0.00020851234 

Nm 

kuxi 

0.00038854941 

Nm 

kMY2 

0.0095461681 

Nm 

k]VLZl 

0.00022903467 

Nm 

kMZ2 

0.00016170646 

Nm 
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Table  C4  Test  point  specific  example  data 


Variable 

Symbol 

Value 

Units 

Gross  strain  gauge  balance  readings 

R 

[0, 

-0.1400, 

0.5400, 

0.1200, 

0.3500, 

0.0800] 

mV/V 

Parent  aircraft  angle-of-attack 

0g,ac 

4.99 

degrees 

Carriage  axes  pitch  offset 

Oac,c 

-3 

degrees 

Carriage  axes  roll  offset 

fac,c 

5.5 

degrees 

Absolute  model  position,  from  store 
support  feedback 

[xg,ss/  ygfss,  Zgfss] 

[6.838,  -109.636, 
175.804] 

mm 

Relative  model  position,  from  store 
support  feedback 

[xrg,ss,  yrg,ss,  zrs,ss ] 

[3.642,  -37.926, 
103.808] 

mm 

Model  orientation  in  tunnel  axes, 
from  store  support  feedback 

\fg,SS,  Og,SS /  y/g,ss\ 

[-5.004,  -8.500, 
0.010] 

degrees 

Tunnel  static  pressure 

Ps 

75.04 

kPa 

Tunnel  total  pressure 

Pt 

126.91 

kPa 

Tunnel  total  temperature 

Tt 

308.502 

degrees 

Kelvin 

[XI]  = 


1  0  0 
0  1.34957955824748E  -  02  1.85872010975616E-04 

0  -  2.01270863980661E  -  04  E33297098851551E -02 
0  -  2.92097108131 194E- 03  1.2276290588619E-03 

0  -  7.31078706779773E  -  04  4.85721 196061 138E  -  03 
0  -  3.78793090574685E  -  03  6.00071467638071E-04 


0 

9.10153188 1 77356E  -  05 
1.411 2649989474 1 E  -  04 
2.49931437883246 
1.061941 11329936E- 03 
3.62245421958098E -  03 


0 

5.5107087643103E  -  05 
3.0659443233 1226E  -  05 
-  9.29033874774701E  -  05 
0.545184191413662 
- 1.568557507596 19E  -  02 


0 


- 1 .33844648392573E  -  04 
- 1.32 1 953389 17485E  -  04 
-1.3228477360321  IE -03 
1.57455876867475E  -  03 
0.546663837059006 
(Cl) 


UNCLASSIFIED 


35 


DSTO-TN-1339 


UNCLASSIFIED 


[D]  = 


1  0  0  0  0  0 
0  1  0  0  0  0 
0  0  1  0  0  0 
0  0  0  1  0  0 
0  0  0  0  1  0 
0  0  0  0  0  1 


C.2.1  Buoyant  Tares 

Buoyant  tares  are  calculated  from  the  averaged  zero  readings  taken  at  the  start  and  end  of 
the  run: 

F XBALT  =  —k  .\  Sln  #0 


F 'ybalt  ~  ky  °o  A 

=  0.000660016068866 

Fzbalt  =  cos  00  cos  to 

=  3.784469738615505 

^  XB ALT  ~  ^ MX 2  COS  ^0  C0S  to  _  ^  MX  i  COS  ®0  to 

=  0.000207024519512 

M  YBALT  —  ~^MY2  00  ~  ^  MY  I  OOS  Oq  COS  (j)(] 

=  -0.000388893127716 

M zbalt  ~  kMZ2  cos  0(]  sin  (j)i}  +  kMZi  sin  0(] 

=  0.000000036784027 


The  equivalent  voltage  for  these  forces  and  moments  is  calculated  and  subtracted  from  the 
gross  strain  gauge  zero  readings  to  give  the  buoyant  zero  readings  (Ro).  Note  that  the 
strain  gauge  balance  calibration  equation  is  simplified  due  to  the  first  order  calibration 
being  used. 
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Rr  =  [D\'[X  l]/f„ 

0 

-0.000753116592538 
0.050444371598445 
0.001051155549775 
-  0.000095972088956 
-0.000501243687842 

R-0  =  ^0  "  R-T 

0 

0.000933117886538 
-0.049674565073445 
-0.000691651292775 
-0.000239929066044 
0.000531765222842 _ 

C.2.2  Gross  Aerodynamic  Loads 

For  each  test  point,  the  strain  gauge  balance  readings  are  corrected  for  buoyant  tares.  The 
gross  aerodynamic  loads  are  then  calculated  from  these  readings: 

R  =  R  -  R0 

0 

-0.140933117886538 
0.589674565073445 
0.120691651292775 
0.350239929066044 
_  0.079468234777158  _ 

Ha=[Xl\'[D\R 

0 

-9.697177517617988 

44.126648133526913 

0.045657007450587 

0.645302229378291 

0.151917836927412 
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C.2.3  Sting-Balance  Deflections 

Deflections  of  the  sting-balance  assembly  are  calculated  using  gross  loads: 

kw  0  0  0  0  0 

0  ks  0  0  0  k9 

0  0  k6  0  k-j  0 

0  0  0  ^00 

0  0  k2  0  £3  0 

0  k4  0  0  0  k5 

0 
0 

1.298906930280108 
0.067065578244168 
-0.352910417100972 
-0.071727255307206 

C.2.4  Store  Model  Position  and  Orientation 

Store  model  position  and  orientation  is  measured  by  the  store  support  system  transducers. 
Positions  and  orientations  are  output  in  tunnel  axes  referenced  both  to  the  tunnel 
coordinate  system  origin  ("absolute")  and  to  a  user  specified  location  ("relative"),  which  in 
this  case  is  the  carriage  position.  Model  orientations  and  both  absolute  and  relative 
positions  are  corrected  for  sting-balance  deflections.  Model  orientations  and  relative 
positions  are  then  transformed  into  carriage  axes,  and  the  radial  distance  and  polar  angle 
are  calculated. 


The  following  Euler  rotational  transformation  matrices  are  required: 

=  m,  (-  ^  ]u,_  (-  U,  (-  J 

cos  0gja  cos  Vg,ss  sin  sin  eg,ss  cos  Vg'Ss  -  cos  tg*  sin 
=  cos  Ggss  sin  y/gss  sin  </>gJU  sin  Ggss  sin  y/g,ss  +  cos  </>gss  cos  i//gss 


iS 

sin^  cos  6 

1  g,ss  g,ss 

cos 

$g,s: 

,  Sln^„, 

,<x>WgA 

s+Sm4,s 

s  Sin^ 

cos 

: _  sin  </>g  ss 

COSl// 

cos  d)  cos# 

'  g,SS  g,ss 

0.9852463  0.0871996  -0.1472613 

-0.0862672  0.9961908  0.0127189 

0.1478094  0.0001726  0.9890158 
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cosrjcosz  sin  k  sin  ^  cos  cos  v  sin  x  cosvsinz/cos  %  +  smvsin% 

=  cos^sin^  sinvsin^sin  j  +  cosvcos^  cos  v  sin  rj  sin  x~  sin  v  cos  % 
-sin  7  sin  v  cos  77  cos  v  cos  7 

"  0.9999802  0.0012447  -0.0061609" 

=  -0.0012519  0.9999985  -0.0011628 
0.0061594  0.0011705  0.9999803 


0.9852463  0.0871996  -0.1472613  0.9999802  0.0012447  -0.0061609 

-0.0862672  0.9961908  0.0127189  -0.0012519  0.9999985  -0.0011628 

0.1478094  0.0001726  0.9890158  0.0061594  0.0011705  0.9999803 


0.9842107  0.0882534  -0.1534297 

-0.0874342  0.9960969  0.0120917 

0.1538980  0.0015142  0.9880856 


Ms,c  =  M,{-es  J 

cos  9gac  0  sin  9gac 
=  010 

S>n  dK,ac  0  COS0g,ac_ 

~  0.9945219  0  0.1045285" 

=  010 
-0.1045285  0  0.9945219 

cos  9ac  c  sin  <pacc  sin  9acc  cos  <pac  c  sin  9ac  c 

0  cos  <j)aC'C  -sin  <j)acc 

sin  °ac,  sin  <j)acc  cos  9acc  cos  <j)acc  cos  9acc 
"0.9986295  -0.0050162  -0.0520950" 

=  0  0.9953962  -0.0958458 

0.0523360  0.0957144  0.9940320 

Mg,c  =  Mg,acMac,c 

'  0.9993969  0.0033283  0.0345652  ' 

=  0  0.9953962  -0.0958458 

-0.0347251  0.0957879  0.9947959 
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Mg,b  =  MgACMacfMc,  =  M  gcMcb 


Mc,b  =  MlMg, 


0.9782730 

-0.0690144 

0.1954968 


0.0881476 

0.9919498 

-0.0909148 


-0.18764861 

0.1061721 

0.9764812 


Store  model  position  is  calculated  by  transforming  strain  gauge  balance  deflections  into 
tunnel  axes  and  adding  them  to  store  support  system  encoder  feedback: 


1 

<50 

H 

1 _ 

-  x'~ 

yg> 

= 

y  g,ss 

+  Mg,ss 

y 

_zg,b. 

_Zg,ss  _ 

z' 

6.646721 

109.619479 

177.088640 


1 

X 

1 _ 

1 

X 

1 _ 

~xr 

= 

yrg,ss 

y 

_zrs,b_ 

1 

V* 

N 

_ 1 

z' 

3.450721 

-37.909479 

105.092640 


Store  model  orientations  in  tunnel  axes  are  found  from  the  elements  of  Mg/b: 

Mg  >fc(3,l)=  -sin  6gh  =0.1538980 


9gb  =  -8.8529° 

Mg,h (3,2)  =  sin  (j)g  b  cos  9gb  =  0.0015142 
M  gb (3,3)  =  cos  (j)g  b  cos  9gb  =  0.9880856 


^g,b 


f 

=  arctan 

v 


MjVf' 


=  0.0878° 
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Mgb(\,\)  =  cosOgb  cosy/gJ)  =  0.9842107 
Mgb( 2,\)  =  cos6gb  sin y/gh  =  -0.0874342 


( 

y/gb  =  arctan 

V 


MjV) 

MjU) 


=  5.0767° 


A 


J 


Likewise,  store  model  orientations  in  carriage  axes  are  found  from  the  elements  of  Mc,b- 

Mcb{3,i)=-sin0cb  =0.1954968 


0cb  =-11.2737° 


Mcb{ 3,2)  =  sin  <j)cb  cos  6cb  =  -0.0909148 
Mc  b  (3,3)  =  cos  (f)c  b  cos0c  b  =  0.9764812 


Kb  =  arCtaI1 


=  -5.3192° 


Mc  b(lX)=cos0c,b cos  Vc,b  =  0.9782730 

Mcb( 2,l)  =  cos  6^  sini//cb  =  -0.0690144 


/ 

y/cb  =  arctan 

v 


mJ»)' 

J 


=  -4.0354° 


Model  positions  in  carriage  axes,  relative  to  carriage  position,  are: 


1 

H 

1 _ 

1 

1 _ 

yrc,b 

yrg,b 

1 

N 

_ l 

1 

N 

_ 1 

-0.200709 

-27.656858 

108.298462 


The  grid  radial  calibre  position  and  lateral  polar  angle  can  then  be  calculated.  Care  should 
be  taken  that  yrc,b,  zrc,b  and  D  all  have  consistent  units.  In  the  following  calculation,  mm  are 
used. 


r/  _ 

/D~  D 
=  4.5492 
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O  p=  arctanl 


\Zrc,b  ) 


=  -14.3258° 


C.2.5  Aerodynamic  Loads  in  Balance  Axes 


The  corrected  store  model  orientations  are  used  to  calculate  wind-on  metric  mass  tare 
loads,  which  are  then  subtracted  from  gross  loads  to  give  aerodynamic  loads  in  balance 
axes: 

^ XBALT  =  —k-X  S'n  @g,b 

=  0 

FyBALT  =  COS  Qg.b  ^g.b 

=  0.00566317 

F ZBALT  =  k-z  COS  @g,b  C0S  ^g,b 

=  3.73938002 

^ XB ALT  =  k.UX 2  COS  ®g,b  C0S  ^g,h  ~  COS  ® g,b  S'n  ^g.h 

=  0.00019326 

kMY i  cos0g,b  cos4 

=  0.00108522 


^ yb alt  ~  kMY2  sin  6gb 


rgfi 


M  zb  alt  ~  k\fz  2  cos  0gb  sin  <j)gb  +  kMZi  sin  0gb 

=  -0.00003500 


ht  = 


XBALT 


YBALT 


F. 


ZBALT 


M 

M 

M 


XBALT 


YBALT 


ZBALT  . 


H  =  Hr-H1 


0 


-9.70284069 

40.38726811 

0.04546375 

0.64421701 

0.15195284 
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C.2.6  Aerodynamic  Loads  in  Body  Axes 


Rotational  misalignment  between  the  balance  and  body  axes  is  represented  by  the  non- 
Euler  transformation  matrix  Mb,bd-  most  cases,  any  small  offset  will  be  nulled  during 
initial  model  setup.  If  there  is  a  large  misalignment,  care  must  be  taken  to  properly 
account  for  the  misalignment  at  various  stages  in  this  process.  This  case  is  beyond  the 
scope  of  this  document.  In  this  example  there  is  no  misalignment,  so  the  transformation 
matrix  is  the  3-dimensional  identity  matrix. 


M 


b,bal 


1 

0 

0 


0 

1 

0 


0 

0 

1 


The  origin  of  the  store  body  axes,  the  Model  Moment  Reference  Centre  (MRC),  is 
displaced  from  the  origin  of  the  balance  axes,  the  Balance  Moment  Centre  (BMC)  by  the 
vector  s mrc, bmc-  The  aerodynamic  loads  about  the  store  body  axes,  F  and  M,  can  be 
calculated  from  the  aerodynamic  loads  about  the  balance  axes,  Fbal  and  Mbal,  as  follows: 


s’” 

II 

II 

K> 

Lj 

co 

uT 

ll 

ll 

s 

F  = 

M  bjyaiF BAL 

0 

= 

-9.70284069 

40.38726811 

~xi~ 

0.0025' 

S  MRP,  BMC  ~ 

I 

_ l 

— 

0 

0 

M  ~  LSbalM Bal  + 


0 

-z. 


yb 


o 


-x. 


-yb 


o 


F 


0.04546375 

0.54324884 

0.12769574 


Note  that  when  translating  the  moments,  the  forces  in  store  body  axes  are  used.  The 
displacement  vector  should  therefore  be  in  store  body  axes.  Any  deflection  between  the 
BMC  and  MRC  is  accommodated  during  balance  calibration,  and  it  is  assumed  that  the 
assembly  is  rigid  between  the  metric  end  of  the  balance  and  the  MRC. 
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C.2.7  Tunnel  Flow  Conditions 
Tunnel  Mach  number  is  calculated  as: 


=  0.90 


Dynamic  pressure  is  calculated  as: 

q  =  —  PM2 
2 

=  42.544  kPa 

Reynolds  number  is  not  used  for  further  results  in  this  example,  and  is  not  calculated  here. 
C.2.8  Aerodynamic  Force  and  Moment  Coefficients 

The  aerodynamic  loads  in  body  axes  are  used  to  calculate  body  force  and  moment 
coefficients. 


Cx 

CY 


Ijl 

qS 

Fj_ 

qS 


0 

-0.4812 


Cz 

C, 


=  F^  =  2.0028 
qS 

=  ^  =  0.0918 
qSd 

=  —L  =  1.0964 
qSd 

=  — =  0.2577 
qSd 


The  body  force  and  moment  coefficients  can  be  transformed  into  missile  axes,  or  non¬ 
rolling  axes,  using  the  Euler  rotational  transformation  matrix: 

"l  0  0 

Mpb  =  0  cos^c>,  -sin^ 

_0  sin  <f>ch  cos<pcb 
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1 
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ii 

Cj 

Cy 

L  zp  J 

Cz_ 

0 

-0.2934 

2.0387 


1 

.O 

_ l 

1 

1 _ 

Cmp 

II 

1 

1 _ 

_c„_ 

0.0918 

1.1156 

0.1550 
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